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Weak gravitational lensing of background galaxies offers an excellent opportunity to study the 
intervening distribution of matter. While much attention to date has focused on the two-point 
function of the cosmic shear, the three-point function, the bispectrum, also contains very useful 
cosmological information. Here, we compute three corrections to the bispectrum which are nominally 
of the same order as the leading term. We show that the corrections are small, so they can be 
ignored when analyzing present surveys. However, they will eventually have to be included for 
accurate parameter estimates from future surveys. 
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O ' I- INTRODUCTION 

(N : 

C3 ' Weak lensing offers cosmologists the opportunity to probe the distribution of mass in the universe Q . This prospect 
is so alluring because theories make first-principles predictions about this distribution, so we can hope to extract 
' important constraints on fundamental cosmological parameters from weak lensing surveys ^3 ■ 

\l In many senses, this promise is similar to that felt by those who studied the cosmic microwave background (CMB) 
^ ' a decade ago: theoretical predictions are straightforward; experiments have detected the effect (anisotropics in that 
, case and cosmic shear in this); and there are grand plans for the future (which have been realized in the case of the 
Cn : CMB). 

[ Armed with this optimism, cosmologists are quick to throw in warning labels: the signal is extremely small, so the 

• only way to measure cosmic shear is to average over many background galaxies. Each individual galaxy is observed 
' with its own set of systematics (seeing, elliptical point spread functions, calibration, unknown or at least uncertain 
. redshift, etc.) and these vary from one galaxy to another. Measurements of cosmic shear are unlikely to produce 

smooth maps since there inevitably will be bright stars which must be masked out. Accounting for these masks leads 
to complicated window functions. It is not clear then that weak lensing measurements will eventually pay off as did 
those of the CMB. 

There is one area though in which weak lensing measurements have an advantage over the CMB: the higher point 
Q functions of the cosmic shear field are potentially simpler to interpret and more relevant than those in the CMB. 
i—{ _ Naively, this is what one would expect, for the cosmic shear field is sensitive to the matter density which has gone 
c/2 . nonlinear and therefore will have large corrections to the Gaussian limit. Temperature fluctuations on the other hand 
' are still stuck at the 10~^ level, so are expected to be very close to Gaussian (recall that the initial distribution of both 
^ . temperature and matter inhomogeneities was likely Gaussian and for a Gaussian distribution the higher point functions 

• are trivially related to the two-point function). We might expect then the 3-point function of CMB anisotropics to be 
rS I very small, while that of the cosmic shear to be quite large, at least on small scales. Much work has been done over 

■ the past few years attempting to debunk these naive ideas. We have found that the higher point functions of the CMB 
are quite interesting: lensing TH IT^ hot gas 0|, peculiar velocities psj, and reionization 0,^3 leave their 
imprint in these higher point functions. Nonetheless, the fact remains that to date there has been no detection of a 
non-zero 3-point function, for example, in the CMB,jWhile several groups 0,^3 have claimed such a detection in the 
cosmic shear field. Further, a number of authors IMS S US 01 have showed that the bispectrum of the cosmic shear 
field will be able to constrain important cosmological parameters, including properties of the dark energy. Certainly, 
then, we need to obtain accurate predictions of the bispectrum of the shear. 

With this in mind, it is important to emphasize just how important "higher order" corrections to the bispectrum 
might be. To understand this, recall that, to lowest order, the shear is a line-of-sight integral over the matter 
overdensity: 



7(0)= j dxfix)Si^ix,0)) (1) 

where 6 is the angular position on the sky, / is a weighting function, and x is the comoving distance along the line 
of sight. Roughly then the bispectrum, which is proportional to (7'^), is proportional to {S^). It vanishes therefore in 
the large-scale limit where the density field is Gaussian and nonlinearities are irrelevant. The main contribution to 
the bispectrum then comes from the fact that, due to gravity, S evolves nonlinear ly. In perturbation theory, we would 
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write 5 = S^^") + S'^^^ + . . . where S^^'' is the hnear overdensity and (5*^^^ is proportional to 5j^. This main contribution 
to the bispectrum then comes from terms proportional to {S''-^^ S^^'^ S^"^^) and so is proportional to Sj^. 

It is clear then that any correction which alters the linear relation between shear and overdensity in Eq. is of 
the same order in Sl as the "main contribution." Here we study three such corrections first identified by Schneider 
et al. and compute their effect on the bispectrum: 



• Reduced Shear We estimate shear by measuring ellipticities of background galaxies, invoking the relation |lLl2,'^| 
e-i = 2ji, where the subscript i refers to the two components of ellipticity/shear. This relation though is only 
approximate; the full relation is 

where k is the convergence. Expanding the denominator, we see that the ellipticities used to estimate cosmic 
shear have terms quadratic in the perturbations (27^^). This quadratic term contributes to the bispectrum at 
the same order as the main term. 

• Lens-Lens Coupling Eq. does not account for the fact that lenses are correlated along the line of sight. 
This lens-lens coupling induces another quadratic term in the relation between shear and overdensity. 

• Born Approximation When computing the shear one integrates along the photon path back towards the 
source. There is a complication inherent in this integration encoded in the argument of the overdensity in 
Eq. (^3): what is the position of the photon at radial distance x if its observed angular position today is 0? 
the naive answer x*^°) — 1) is the position corresponding to the path taken by an undeflected photon. 
Expanding about this zero order position leads to a correction proportional to V(5 ■ [x-x(o)]. Since X differs 
from the undeflected position only if 6 is nonzero, this correction is second order in 6. It too contributes a term 
of order 5'^ to the bispectrum. 

The next section reviews some basic lensing results including the standard computations of the power spectrum 
and the bispectrum. §111 computes the correction to the bispectrum from the three effects enumerated above. The 
bispectrum cannot be simply plotted on a 2D graph since it depends on the three variables required to specify a 
triangle. Therefore, §IV examines various ways of condensing the information contained in these corrections. The 
goal is to see whether the corrections are important. 

II. REVIEW OF BASIC RESULTS 

The deformation tensor is defined as the deviation from unity of the Jacobian relating the undeflected position (^s) 
to the actual position {0): 

V-^i = % - -g^- (3) 

The elements of this 2x2 matrix are the two components of shear and the convergence: 

These are simply definitions. The physics comes from solving the geodesic equation and expressing the distortion 
tensor in terms of the gravitational potential $: 

f^' d'^^(x(0 y)-y) 

V'., {0, Xs) = dxWix, Xs) gJ-^g^J [Skj - i'k, {0. X)] • (5) 

Here, we are assuming a flat universe; sources are assumed to lie at Xs] the first argument of the $ is the 3D comoving 
position (in the small angle limit) 

x(0, x) = [xe.^xOy, x] - / dx'W{x', x)4 v<i> (x(^?, x'); x) (6) 

Ja X 

while the second (x) refers to the cosmic time at which the photon path passed by this position. The weighting 
function is 

W{x. Xs) = 2x(l - x/xs)e(xs - x). (7) 
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We observe ellipticities of background galaxies e^, which are related to the elements of the distortion tensor via 

Q = -j ■ 

1 — K 

One way to estimate the convergence, which is the projected density, is to work in Fourier space; then, 



(8) 



= ^r,(l)e~,(l). (9) 

Here the variable conjugate to 9 is 1. As usual, large I corresponds to small scales. The trigonometric functions in 
Eq. © are defined with respect to an arbitrary x- axis as 

Ti(l) ^ i^ = ^cos(20,) 
2 

Here is the angle between 1 and the x- axis. In the limit in which 2^i, the estimator in Eq. @ reduces to 



T2(\) = IJy = -sin(20i). (10) 



reduced shear 

^(1) cos(2,/.;)7i(l)+sin(2,/.;)72(l) 

cos(20O(Viii(l) - fe(l))/2 + sin(20O^i2(l) 

lens— lens ,2 

^ —m- (11) 

The approximation on the first line neglects the fact that ellipticities are sensitive to the reduced shear; the approx- 
imation on the third line neglects the second order term in the brackets in Eq. (Q, the term which accounts for the 
fact that the lens distribution is correlated. The projected potential in Eq. Hll|l is defined as 

^{6, xs) = r dx-^^%^i> (x(^, x);x). (12) 

Jo X 

In the Born approximation the gravitational potential in the integral along the line of sight is evaluated at the 
unperturbed path. In that case, its Fourier transform reduces to 

/-Xs T/f/(' 

V;(1,X.)'^^(°)(1,X.)- / rfx 2 > ftx) (13) 

Jo X 



with 



1 f dks^ 



0(1, X) = ^J k,;x)e'''''. (14) 

Note that (p is dimensionless unlike $ which has dimensions of (length)^. 

The statistics of ^'^^^ follow directly from those of 0, which are simple if we include only modes with small, i.e. 
the Limber approximation ,23. .24] . 

m,xW,x')) = {2n)H^\ + mx^x')P'f{i/x;x)/x^. (15) 

Here P,s> is the power spectrum of the gravitational potential. Similarly, the three point-function is related to the 
spatial bispectrum 16]: 

(<^(li, Xi)<^(l2, X2)<^(l3, X3)) = (27r)2^2(i^ ^ 1^ ^ _ ^y^^ _ x")B^{h/x, h/x, h/x; x)/x'- (16) 

Then, we have 

(^W(l)V;(°Hl')) - {27rfS^{l + l')P2il), (17) 

with 

P2(0 = y^ dx ^^, ' p^{i/x;x)- (18) 



Similarly, the three-point function is 

(Vi(°)(li)^(°)(l2)^(°)(l3)) = (27r)2^2(i^ + 12 + 13)^3(1^,12,13), 
where now the projected power is a line-of-sight integral over the bispectrum: 



P3(ii, 12,13)= r 

Jo 



dx- 
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-B-s-ih/Xih/Xih/x^x)- 



(19) 



(20) 



The power spectrum, or the Cfs, are defined as the coefhcient of (27r)^5(l + 1') when computing the variance of k. 
Since k = —Pij^^^/2 in the standard computation, we have 



Ct = /^P2(/)/4. 

Similarly, the bispectrum of the k estimator is 

B«(ll,l2,l3) = -/^P3(ll,l2,l3)/8, 



(21) 



(22) 



in agreement^ with previous results The bispectrum with all /'s equal, the equilateral configuration, is shown in 
Fig.d 
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FIG. 1: Equilateral bispectrum. Solid curve is the standard result; short-dashed [red] curve is order of magnitude of corrections 
considered here; long-dashed [blue] curve is cosmic-variance error. The signal to noise from a single configuration is therefore 
extremely small. 



^ One subtlety when comparing with other results is the sign. The sign here is negative because # oc —5. 
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A few qualitative comments are in order here. The standard measure of the ampUtude of fluctuations is Z^C" = 
Z^P2/4. Let's do an order of magnitude estimate for this quantity l^P2 in terms of the ampHtude of density fluctuations, 
A2 = k^Ps /2Tr'^. Since I ~ (k/Ho) and since $ ~ {Ho/kf5, we have P$ ~ Ps/l'^- Now, Eq. ^ suggests that 
P2 ~ P<s>/x^\ since x ~ l/k, P2 ~ {k/lfP^ - {k/lfPs/l'^ - AVZ^. We expect then that l'^P2 should be of order 
A^/L What this means physically is that projection effects suppress the 2D power spectrum by a factor of l/l. There 
is even a nice explanation of this in terms of Fourier modes [IS 113: only modes with small fca contribute; these are 
a fraction of l/l of the total number of modes. The bottom line then is that the power spectrum of the convergence 
is smaller than the power spectrum of the 3D density field. 

Similar order of magnitude estimates relate the angular bispectrum, —l^P^/S, to the 3D bispectrum of the density 
field: 

l^Ps ~ ^ H^Bs. (23) 

The corrections we consider below are all of order PP-^, which by the arguments of the preceding paragraph are 
of order HqP^. The 3D bispectrum Bfj is nominally of the same order as the square of the power spectrum, P|. 
However, numerically it is a bit larger |20|, as indicated in Fig.^ so the corrections we compute are not as important 
as we would have hoped. 



III. HIGHER ORDER TERMS 



We now compute the corrections to the bispectrum from going beyond the approximations in Eq. and Eq. H13|l . 



A. Reduced Shear 



The first-order correction to the e — 7 relation is 



(24) 



When we switch to Fourier space, the relation between ellipticity and shear is a convolution integral (products in real 
space correspond to convolutions in Fourier space): 



^T(i) 



(2^ 



■r2r,(i-i')^(")(i')^(°)(i-i') 



(25) 



When we form the bispectrum of the n estimator (Eq. ©), the first order correction emerges by replacing one of the 
three ellipticities with the higher order Eq. H25|) . So this correction to the bispectrum estimator becomes: 



(^(Il)^(l2)^(l3))^ 



Alf J {2n 
(li ^ I2) + (li ^ I3 



rl'^Tdh - l')(^(°'(l')^'"'(li - l')i'^"Hh)i'^"\h)) 



(26) 



The four-point function for the potential gets contributions from the connected part - the trispectrum - and the 
disconnected part: the product of power spectra. Here we consider only the latter set of terms as these are expected 
to dominate. That is, let 

X "(52(1' + 12)^2(1^ _ J/ + 1,^) + ^2(j/ ^ j^)^2(i^ j^^ 



The integral over 1' then leaves the coefficient of (271)^5^(11 + 12 + I3) as 

mi)qP2{i2)iiP2{i3) 



4/2 



[qT,ih) + iim2)] 



plus permutations. The contraction over the geometric factors involves 

T.(li)r,(l2) = Mcos(20i2) 



(27) 



(28) 



(29) 
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where (j)i2 is the angle between li and I2. 

Defining -S123 = ^(li, I3) as the coefBcient multiplying (27r)^(5^(li + I2 + I3), we therefore have 

i?r23 = ^^'^^'^'^jg^^'^^'^ [COS(2013) + COS(2012)] (30) 

plus permutations. 



B. Lens-Lens Coupling 



Lens- lens coupling in encapsulated by the second term in square brackets in Eq. Q. This contribution to the 
distortion tensor is then 

^lb{e, xs)^-j dxwix, xs)'^,acAb{e, x)- (3i) 

At second order in this reduces to 

^Ue. Xs)^- I dxWix, Xs)'f,ac{d, x) I dx'wix', x)'fMO, x')- (32) 



Using Eq. @, we can compute how this term contributes to the estimator of convergence. The second-order 
contribution is 

AM) = 

m) f dH' f ^w{x,Xs) I ^M/(x',x)<^(l',x)0(l-l',x')- (33) 



12 J (2,^)2 ' ' J ^2 ^,2 

Here the geometrical factors are defined as 

-El(ll,l2) = ll ■ I2 [ll,xl2,x — ll,j/l2,y] 

i;2(ll,l2) = \l-h[Kxh,v+KyhA- (34) 

The estimator for the bispectrum of the convergence then gets a Gaussian contribution from the lens-lens term. 
One such term is 



{k\\,)ki^Hh)k^^\x,)) = zMM^ / %w{x,xs) f %wix',x) I ^wix2,xs) f %mx3,x.) 

X J (02 ^^(1'' li - l')(^(l', X)^(li - 1', X')0(l2, X2)0(l3, X3)). (35) 

Two other terms exist with li ^ I2 and li ^ I3. In this case, the trispectrum does not contribute in the Lim- 
ber approximation. Physically, the Limber approximation sets all lenses close to each other; mathematically, this 
corresponds to enforcing the constraint that the line of sight distances are all equal Here this constraint sets 

X = x' = X2 = X3i so that W(x',x) vanishes. The only relevant terms therefore are the two products of two-point 
functions. Momentum conservation from the first such pair enforces 1' = — 12 and li — 1' — 13. Each two point 
function is evaluated using the Limber approximation as in Eq. p5|) . 
The bispectrum from lens-lens coupling is then 

-S123 = j^2 y (x,x.) y ^w^(x ,x)w^(x ,xs) 

X [P$(/2/x; x)^*('3/x'; x') + i^*(^3/x; x)^'*('2/x'; x')] + (ii ^ I2) + (h ^ h)- (36) 

Here we have used the fact that i?i(li, — 12) — £'i(li, I2) — i?(l2, li). 
The geometrical factor in front reduces to 

T.(10£;.(l2, 13) = ^'^'^'7^'^^^^ cos(<^i2 + ^13). (37) 
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Therefore, 



X" J X" 

X [P$(/2/x; x)^*a3/x'; x') + i^^Cs/x; x)^*a2/x'; x')] + (ii ^ I2) + (h ^ I3). (38) 

C. Born Approximation 

The distortion tensor in Eq. © evaluates the potential everywhere along the unperturbed path of the light. To go 
beyond the Born approximation, we need to evaluate the potential at x = Xq + i5x where 

Sxa{e,x) = - f dx'w{x\x)^'^>A^o;x')- (39) 

This leads to a new contribution to the distortion tensor, which in Fourier space, reads 

^^'^(i, X.) = - / ^wix, xs) I ^wix', x) / -^O'tUi - i')c0(i'; xm - 1'; x')- (40) 

This extra term in the distortion tensor contributes to the estimator for the convergence 

.Born ^ -m J ^w{x,x.) J ^W{x',x) J ^F,(iM-i')0(i';x)^(i-i';x'). (41) 

Here 

F,(li,l2)=21i -127,(11). (42) 

This expression is identical in form to that from lens- lens coupling, with the substitution Ei — > ~Fi. We can therefore 
copy the result from Eq. to get 



Bf- - / ^H^^(X,X.) / ^I4^(X',X)T^(X',X.) 

X [P*(;2/x;x)i^*(^3/x';x')^^Ki2,i3) + P*(?3/x;x)^^*(^2/x';x')^.(i3,i2)] + (h ^ I2) + ih ^ h)- (43) 

But, r,(li)F, (12,13) = Z?ZiZ3Cos(20i2)cos(</)23), so 

^Bo.„ ^ I !^w^u,x.) I '^wix\x)Wix',xs) 

X [Z2'^*('2/x; x)i^*(^3/x'; x') cos(20i2) + ?3'^*('3/x; x)^^*('2/x'; x') cos(20i3)] + [h ^ h) + (ii ^ h)m 

D. Summary 

Here we collect the results from the previous three subsection. The reduced shear correction can be expressed in 
terms of the 2-point function P2 : 

5. ^ MMM^[eos(20i3)+cos(20i2)] 
lb 

'l'3-f2('l)^2(^3) r , N , /„ , XI 

fcos(2023) +C0S(2(/)12)] 



16 

liltP2{l2)P2{h) 

16 

The other two corrections are best expressed in terms of 



[COS(2013)+COS(2023)] (45) 



i{hM) ^ m / ^w^'(x,x.) / ^mx',x)mx',x.)^*(^i/x;x)i^*(^2/x';x')- (46) 
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Then, the lens-lens term is 



Rll -hh COs(2(/)23) C0S(.^12 + 013) r^^, , n , w, ,M 
^^123 = ^ [J(,'2,i3) + J(,t3,i2)J 

;1;3COS(2013)COS((/)21 + fe) rw, , N , w, ,M 
- ^ [-'(,'l,t3j +-'l'3,tljj 

_ /2^1COS(202O COS(032 + 031) j,^^^^ ^ ^^^j ^ ^^.^ 



And the Born term is 



^Born ^ '^""^('^23) ^^2^^^^^ C0s(2</)12) + /2/(Z3, ^2) COs(20i3)] 

- H^^^ ^3) COS(2012) + lll{h, h) COS(2023)] 

- H^^^ [/2'/a2, ^l) COs(2</)32) + ^2) COS(2013)] . 



(48) 



One word of caution: all of the above assume that our estimator for k, is as given in Eq. 0. One could also 
imagine defining the bispectrum as the three-point function of one-half of the trace of the distortion tensor. These 
two expressions agree in the zeroth order case, but they disagree when these higher order corrections are included, 
because the distortion tensor is no longer the second derivative of a potential ip. Were we to be interested in the latter 
definition, then the cos(2(/)) terms inside the square brackets in Eq. (|48|l would be replaced by 1. For the lens-lens 
term, the product of cosines on the first line of Eq. (|T7l) would be replaced by cos^ 023 ; on the second line by cos^ 0i3 ; 
and on the third by cos^0i2. Practically, we think that the estimator of Eq. ||5J) is more relevant, since it is the 
different components of ellipticity that are measured, not the distortion tensor. 



IV. HOW IMPORTANT ARE THE CORRECTIONS? 



There are a number of ways of assessing the importance of the corrections considered here. First, we compute the 
skewness as a function of smoothing angle. The smoothness at any given angle is an integral over all configurations of 
the bispectrum with a particular weighting scheme. Thus it reduces all elements to a single number. Second, we can 
compute the equilateral configuration as a function of multipole moment I. Comparing the signal with the anticipated 
noise allows us to see whether this one configuration, in which all Z's are equal, is sensitive to the corrections. Finally, 
we compute the anticipated bias on a cosmological parameter from future measurements of the bispectrum if these 
corrections are neglected. If this bias is very small, smaller than the anticipated statistical error, then there is no need 
to worry about the corrections. 



A. Skewness 



The convergence skewness is defined as 

^3 ^ {^)/{^? (49) 

where R is the convergence smoothed over certain window function. In the weakly nonlinear regime where second 
order perturbation theory applies, Bernardeau et al. 27] showed that the skewness does not depend on the density 
fluctuation amplitude cg but is very sensitive to the mean matter density This behavior holds even in the highly 
nonlinear regime ,2\\ . So 6*3 is particularly useful to break the degeneracy of f2„i and erg in the lensing power spectrum. 
Detections of skewness have been reported by several groups 0, ; future surveys such as Canada- France-Hawaii 
Telescope Legacy Survey "2^ could determine ^Im to 10% in this fashion. 

Corrections to the lensing bispectrum affect the prediction of skewness and thus bias the constraints of cosmological 
parameters. We quantify corrections of reduced shear, lens-lens coupling and deviation from Born approximation to 
{k^). {R?) is related to the lensing bispectrum by 

{r') = J B{h,h,h)Wih)W{h)W{k)^-^ . (50) 
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Here, W(l) is the Fourier transform of the window function W{9). We study two window functions, the compensated 
Gaussian and the aperture: 



WcciO) = (1 - 0V2^/) expMV2e/) 



(51) 



where 9f is the characteristic scale in both cases. For both, (k'^) is then a function of 9f. The corrections from the 
three effects considered in §111 and the total are shown in Fig. |2 on interesting scales, corrections are smaUer than 
about 2%. Since 5*3 oc where a ^ 0.8 113, H^, these corrections could bias the determination of flm by less than 
about 2%. So they can be safely neglected in the near future. 
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FIG. 2: Corrections to lensing skewness. The window functions we adopt are compensated Gaussian (top panel) and aperture 
function (bottom panel). Individual and combined corrections are less than about 2% on all scales of interest. 
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B. Equilateral Configuration 

One configuration which is often used as a standard is the equilateral configuration, wherein li = I2 = I3 ^ I. The 
bispectrum can then be plotted as a function of I. Let's consider the corrections to the equilateral bispectrum. 
For the reduced shear correction, all the cosines in Eq. H30|l are —0.5, so adding up all the permutations leads to 

BlI,^^. (52, 

For the lens-lens correction, in the first line of Eq. (|47|l (f>i2 + (pia — 2tt; cos(2(/)23) = —1/2; and all three permutations 
contribute equally, so 

Bfu = —I{l,l). (53) 
In the equilateral case of Eq. H48|) . all the cosines are equal to —1/2, so 

BEr = ~BI\J2. (54) 

The resulting corrections are shown in Fig. |3 

The new terms are only about ten percent of the first order term usually considered in this equilateral configuration. 
Nonetheless, they may still be important for precision cosmology where we sum over many different configurations. 



C. Cosmological Parameter Bias 



There is a simple formula relating the error in a cosmological parameter to a mis-estimate in the theoretical 
prediction. 



Ap = F 



W{li,l2,h) 



dp 



-AB, 



h,l2,h- 



(55) 



Here the bias in the parameter is Ap; F is the Fisher matrix (here just one number since we treat the simple case 
of only one parameter); w is the weight, or the inverse variance, from the experiment of interest, and AB is the 
mis-estimate in the bispectrum, here taken to be the full set of corrections computed above. The Fisher matrix too 
depends on the survey. It is 



(dBi, 



dp 



(56) 



The weights for a particular configuration depend on the survey in question. We are interested in the question of 
whether these corrections can ever be important, so we take the minimum possible errors: cosmic variance due to 
simple Gaussian fluctuations. Following Takada and Jain 4] , the weights are 



^(^1,^2,^3)"^ = Ai23Ci,Ci,Ci, 



(57) 



where the Ci's are l'^P2{l)/4:, and A = 1 when all Vs are different, A = 2 when two Vs are equal, and 6 when all Vs 
are equal. Under this weighting, the sum extends over h < h ^ h- 

A simple application of this formula is to consider the parameter to be the amplitude A of the bispectrum assuming 
the shape is known. Then, Eq. (|55|l reduces to 



/AA 



'^hM,h w{li,l2,l3)Bi^i,j^ABi^ 



bias Eii,i2,i3 ^('l'^2,?3)S;^;^;^ 

This is to be compared with the fractional statistical error. 



(58) 



AA\ 



statistical 



^ wih,l2,k)Bl^i, 



-1/2 



(59) 
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FIG. 3: The bispectrum of the convergence estimator for the equilateral triangle configuration: h — I2 — I3 = I. Thin dashed 
line represents old calculation, which neglected second order terms; thick line includes all relevant terms. Red curve is term 
induced by second order relation between ellipticity and shear: e ~ 2'yK,. Blue curve is the sum of the contributions of lens-lens 
coupling and the Born correction. Note that all second order corrections are negative. Top panel shows fractional difference. 



Before evaluating these sums, we can estimate them. The bias is of order /S.B/B which, when summed over many 
modes, led to percent level changes in skewness. Although the weighting scheme is different here, we might still expect 
/S.A/A ~ 1%. The statistical error is of order SBcos. vm./B for each configuration where both the cosmic variance 
error and the bispectrum are plotted in^ The ratio is seen to be about 10"^. This is reduced by the square root of all 
configurations of order [^^ax]^^^ ^ l*-*^. Thus we expect the statistical error to be of order lO^''. Evaluating, we find 



fAA\ 



= 1.08 X IQ-^ 



jias 



= 1.22 X 10"^ (60) 



statistical 



— 1/2 

for an all-sky survey, in agreement with our estimates. The statistical error scales as f^^y , so we expect it to be 
smaller than the bias for surveys that cover areas larger than /sky = 0.01, or 400 square degrees. 

So at least in principle, these corrections will eventually have to be included. Ignoring them would induce a bias to 
the cosmological parameters up to ten times larger than the anticipated statistical error. Presently, of course, we are 
nowhere near these limits, so we can safely neglect the corrections considered here when analyzing lensing surveys. 
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V. CONCLUSIONS 



We have computed corrections to the bispectrum due to: reduced shear, lens-lens coupling, and the Born correction. 

These corrections arc smaller than the canonical term; this stems from the fact that the spatial bispectrum is larger 
than the square of the power spectrum. While the corrections are small and can be neglected in present surveys, 
when areas as large as 400 square degrees come online, cosmological parameters extracted from the bispectrum will 
be mis-estimated unless the corrections are included. 

This work is supported by the DOE and by NASA grant NAG5-10842. 
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